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We explore, by means of inelastic neutron scattering, the extent of changes in dynamic correlations 
induced by thermal poling of vitreous Si02- The measured vibrational density of states shows an 
excess of modes in certain frequency regions as well as a narrowing of the 100 meV peak. These 
findings indicate that such alterations cannot be ascribed to the appearance of new well defined 
vibrational modes, such as those coming from localized topological defects, but rather arises from 
an increase in ordering in the material as attested in a reduced spread of the inter-tetrahedral 
angles. 



PACS numbers: 42.70,42.65.Ky,63.50.x,78.20.e 

The ability of optical glass fibers to generate visible 
(green) light via Second Harmonic Generation (SHG) of 
intense infrared pump light remained unexplained for a 
long time since the first experimental detection of the 
phenomenon in 1985 El. Two main reasons made such 
phenomenon unexpected. Firstly, a fundamental princi- 
ple derived on symmetry grounds for the second-order 
susceptibility tensor |^| , forbids SHG in an isotropic ma- 
terial, such as amorphous silica. Secondly, a well de- 
fined phase matching between the interacting waves is 
required to allow constructive interference leading to ef- 
ficient SHG. Even if some symmetry breaking mecha- 
nism is assumed, the way in which phase matching waves 
was automatically forced appeared as almost a "miracle" . 
The puzzle was essentially solved five years later by Di- 
anov et al. |jj by invoking the emergence of a spatially 
modulated local dc field, Eq, which, via third-order non- 
linearity (x^)i induces a spatially modulated second or- 
der nonlinearity (x^ °c X^-^o)- This so called Electric 
Field Induced Second Harmonic Generation (EFISHG), 
which is responsible for the frequency doubling, relies on 
a third order nonlinearity which is finite in structurally 
isotropic materials. The dc field is induced by photocur- 
rents originated from ionization of color centers associ- 
ated with impurities in the glass structure. In the sim- 
plest model p], quantum interference between two and 
one photon ionization leads to a modulation of the cur- 
rent, resulting in a phase matching (or more properly, 
quasi-phase matching). Under such conditions the sym- 
metry is broken by the local dc field and no significant 
changes in the structure are needed. 



Given the obvious technological interest of the phe- 
nomenon, a substantial effort has been focused towards 
the achievement of a permanent x and to improve the 
SHG efficiency in glass. Nowadays, the so called thermal 
poling technique, which consists in applying high voltage 
(~ 5 kV) at elevated temperatures (~ 270°C - 300°C) is 
able to provide figures for SHG conversion comparable 
to those shown by some inorganic crystals 0]. It is also 
possible to implement quasi-phase matching efficiently in 
periodically thermally poled fibers || . 

Contrary to what could be expected, there is a lack of 
studies carried out by direct means on structural and/or 
dynamical alterations associated with poling. Most are 
either macroscopic || or mesoscopic H in nature, ren- 
dering scarce information of the underlying microscopic 
physics. Up to now, probably the most direct inves- 
tigation of the phenomenon comes from Raman spec- 
troscopy studies ||,|9|,|l0|] . In fact, significant changes be- 
tween treated (photoinduced) and untreated fibers were 
observed at rather well defined frequencies and associ- 
ated with changes in the density of topological "defects" 
in the Si02 glass structure. Such structural modifica- 
tions were assigned as the main microscopic symmetry- 
breaking mechanism leading to SHG f|[ll]]. A direct in- 
terpretation of the Raman data is, however, hampered by 
the need of previous knowledge of the possible structural 
alterations in order to properly assess the concomitant 
variations in the intricate couplings between the optical 
waves and electronic clouds governing the signal inten- 
sities. As a matter of fact, the aforementioned defects 
model cannot completely explain all the experimental re- 
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suits ||. In spite of this, however, it is sometimes pre- 
sented in the literature as the way of explaining the phe- 
nomenon. 

Our aim here is to investigate the microscopic mech- 
anisms leading to the appearance of the second-order 
non-linear optical response in poled glass, using for the 
purpose a beam of thermal neutrons as a probe which 
directly couples to the nuclei. We thus carried out a 
set of inelastic neutron scattering (INS) measurements of 
the generalized frequency spectrum (also referred as vi- 
brational density of states) covering the frequency range 
where previous reports from optical spectroscopies found 
measurable differences between untreated and poled sam- 
ples. 

We have chosen to work with thermally poled samples 
(nowadays probably the most promising poling technic in 
regarding possible technological applications). The sam- 
ples used in the experiments were made of electrically- 
fused quartz (Infrasil from Hcraeus) with dimensions 
40 x 40 x 0.1 mm. The thickness (~ lOO^m) was cho- 
sen to be the smallest compatible with the polishing in 
order to increase the ratio poled volume versus unpoled 
volume. In fact, after thermal poling pi, the nonlinear- 
ity is located in a region ~ 5/im deep under the anodic 
surface independently of the sample thickness, so that it 
is evident the convenience of using thin samples. A ratio 
of poled/unpoled material of ~ 5% was estimated. The 
thermal poling was carried out at 275°C with an applied 
dc voltage of ~ 5 kV for 40 minutes. The poled sample 
consisted in fifteen plates stacked together. For compar- 
ison a similar stack of unpoled samples was used. 

The INS experiments were carried out on the MARI 
chopper spectrometer at the ISIS pulsed neutron source 
at the Rutherford Appleton Laboratory, Oxfordshire 
(UK). The spectrometer operates in direct geometry 
mode, that is the wavelength (energy) of the incident 
neutron pulse is monochromatized by a set of three ro- 
tating devices (choppers) which select neutrons having 
the required velocity. The energy and momentum anal- 
ysis of the outgoing neutrons, after scattering from the 
sample, is carried out by means of time-of-flight measure- 
ments once neutrons are detected by a large array of 3 He 
counters which span a large angular range. Two differ- 
ent incident energies were used (220 meV and 150 meV) 
which enabled the coverage of kinematic ranges with dif- 
ferent resolution in energy-transfers (the achieved resolu- 
tion varies as a smooth function of AE/E, being typically 
of a few per cent). In all the experiments the samples 
were kept at 20 K in order to reduce the multiple-phonon 
contribution to the spectrum, a quantity which becomes 
substantial in measurements of this kind performed at 
room temperature. The quantity which is directly acces- 
sible to a neutron scattering experiment is the double- 
differential (d 2 er/dSld-E) scattering cross-section given in 
units of barns meV -1 Sterad . For a macroscopically 
isotropic material such as a glass, the quantity of utmost 



importance to specify the dynamics is the Z(oj) spectral 
frequency distribution (or vibrational density of states). 
For samples where the dominant neutron-nucleus scatter- 
ing is of incoherent nature, derivation of Z(u>) from the 
measured cross-section only involves a few corrections. In 
the case of silica, where scattering from Si and O nuclei 
gives rise to coherent (interference) effects, the derivation 
of Z(ui) from the measured data becomes more elabo- 
rated. A number of schemes for averaging-out the coher- 
ent effects have been developed and are known to lead to 
reliable estimates of the spectrum. In our particular case 
we have adopted a procedure described in detail in Ref. 
|]l2| which is based on a phonon expansion of the dynamic 
structure factor S(Q, E) cx (Pa/dfldE and consists in an 
iterative scheme designed to correct a number of instru- 
mental and sample-dependent (multiphonon) effects. On 
the other hand, having a coherently-scattering sample 
is here an additional bonus rather than a nuisance. In 
fact, the single-differential scattering cross-section dc/dfi 
derived by integration over energy-transfers of the mea- 
sured response is directly proportional to the static struc- 
ture factor S(Q) and therefore some hints of structural 
modifications could possibly be derived from the same 
set of experimental data. Besides, the wavevector depen- 
dence of S(Q,E) at a given energy transfer provides di- 
rect information about the phase relationships of the mo- 
tion of specific atomic pairs |r| , since the Fourier trans- 
form of S(Q, E = ctant.) is a dynamical pair correlation 
function Da(r, E = ctant.) which serves to identify the 
atoms taking part in motions at such frequency as well 
as their relative phases. 

Figure |l| displays a comparison between spectral dis- 
tributions for the virgin (untreated) and poled samples 
up to a frequency of 120 meV (w 960 cm -1 ). As stated 
above, such a range comprises the one where the Ra- 
man studies reported changes induced by poling. A mea- 
surable difference in Z(io) between the poled and un- 
poled glasses is readily seen for frequencies about 53 meV 
(« 424 cm -1 ), a band comprising 60 meV (» 480 cm -1 ) 
< E < 95 meV (» 760 cm" 1 ) and the peak at 100 meV 
(rs 800 cm -1 ). Notice that only 5% of the sample was re- 
ally poled and therefore the changes are correspondingly 
scaled down. As can be seen from the lower frame of 
the graph, the differences are outside error bars and this 
supports the statistical significance of such findings. The 
measurements carried out using a larger incident energy 
(220 meV) confirmed the above mentioned changes and 
enabled the exploration of the spectral region extend- 
ing up to 180 meV (« 1440 cm -1 ), which comprises the 
strong double-peak structure [ [L3| with maxima at about 
135 meV (« 1080 cm -1 ) and 150 meV (« 1200 cm -1 ). 
Feeble changes were found within that range of frequen- 
cies. However, the lower resolution and statistics pre- 
vents to discriminate a reasonably clear difference and 
therefore, our ensuing discussion will only concern alter- 
ations below 110 meV. 
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measured optical second order nonlinearity. 

In summary, inelastic neutron scattering is capable 
of monitoring the alterations occurring in fused silica 
at a microscopic level after the material is submitted 
to thermal poling. From the present data (the vibra- 
tional density of states) two main conclusions may be 
inferred. Firstly, the modifications in a frequency range 
encompassing the so called Di and D2 Raman lines are 
consistent with dynamical alterations involving rather 
large groups of atoms. This finding is in contrast with 
the highly localized features expected from a creation of 
three-fold and four-fold ring-like "defects" sometimes as- 
sumed as possible cause of the second order nonlinear 
optical response. Secondly, alterations around 53 meV 
and particularly at 100 meV strongly suggest a higher or- 
der in the treated material implying a reduced spread in 
the inter-tetrahedral angle distribution. Further studies 
are needed to clarify whether these changes can justify 
the value of second-order nonlinearities measured after 
poling. 
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